Astrophysics and Space Science 
DOI 10.1007/s»»»».-«.»-»».»-» 



On the orbital and physical parameters of the HDE 
226868/Cygnus X-1 binary system 
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Abstract In this paper we explore the consequences of 
the recent determination of the mass m = (8.7±O.8)M0 
of Cygnus X-1, obtained from the Quasi-Periodic Oscil- 
lation (QPO)-photon index correlation scaling, on the 
orbital and physical properties of the binary system 
HDE 226868/Cygnus X-1. By using such a result and 
the latest spectroscopic optical data of the HDE 226868 
supcrgiant star we get M = (24 ± 5)Mq for its mass. 
It turns out that deviations from the third Kepler law 
significant at more than 1-sigma level would occur if the 
inclination i of the system's orbital plane to the plane of 
the sky falls outside the range w 41-f-56 deg: such devi- 
ations cannot be due to the first post-Newtonian (IPN) 
correction to the orbital period because of its smallness; 
interpreted in the framework of the Newtonian theory 
of gravitation as due to the stellar quadrupole mass 
moment Q, they are unphysical because Q would take 
unreasonably large values. By conservatively assuming 
that the third Kepler law is an adequate model for the 
orbital period we obtain i = (48 ± 7) deg which yields 
for the relative semimajor axis a = (42 ± 9)Rq (w 0.2 
AU). 

Keywords black holes: individual (Cyg X-1)— stars: 
individual (HDE 226868)— X-rays: binaries— X-rays: 
individual (Cyg X-1) 



intensively studied X-ray source s. After the d etec 



tion of its rapid X-ray variability (jOda et al.lll97ll ) and 
the identification of its optical counter part with the 
09.7 lab supergiant sta r HDE 226868 (jBoltonl [l972 : 
Webster fc MurdinI [l973) . it was considered as one of 
the most likely black hole candidates. The X-ray emis- 
sion in Cygnus X-1 is powered mainly by accr etion from 
the s trong stellar wind from HDE 226868 (jPetterson 
1978h . 

In this paper wc investigate some features of the 
orbital geometry of the HDE 226868/Cygnus X-1 bi- 
nary system in order to check the compatibility of 
the so far obtaine d spectroscopic data of HDE 226868 
( Gies et al.l 120031) with the tight constraints on the 
blac k hole mass m = (8.7zt0.8 ) Mp, r ecently determined 



bv lShaposhnikov fc TitarchukI (120071) who exploited the 



existing correlation (jShaposhnikov fc TitarchukI 120061 ) 
between the low frequency Quasi-Periodic Oscillations 
(LF QPO) and the photon index of the power law spec- 
tral component of the Cygnus X-1 X-ray emission. We 
will look for deviations from the third Kepler law, which 
will not be assumed a priori valid, by discussing the 
conditions at which they occur and their physical plau- 
sibility. 



1 Introduction 

Cygnus X-1, one of the brightest high-energy sources 
in the sky, with an average 1 — 200 keV energy flux 



of 



3 X 10 erg cm s , was discovered by 



Giacconi et al.l (|l967n and soon became one of the most 
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2 Deviations from the third Kepler law 

The analysis of 115 optical spectra of the Ha emission 
fine of HDE 226868 for the 1998-2002 interval allowed 



Gies et al.l (|2003[ ) to phenomenologically determine the 
orbital period P and the projected velocity semiampli- 
tude K (see Table [I|: for P t he same value, within th e 
errors, of the one obtained by Brocksopp et al. ( 1999f ). 
who used a 26-yr data set, was obtained. A pre vious 



19821) 



orbital fit can be found, e.g., in (jCies et al 

From P and K, in turn, it is possible to obtain 
the stellar projected baryccntric semimajor axis xm = 



2 



Table 1 Relevant orbital parameters of the HDE 
22686 8 /Cysnus X-1 binary system. P is the orbital pe- 
riod (jBrocksopp et al]ll999ll. K is the projected velocity 



semiamplitude jcies et alT 20031), q = m/M is t he ratio of 



the black hole to the star mass (iGies et al.ll2003l). m is the 
black hole's mass ijShaposhnikov fc Titarchukll200'ii l. All of 
them have been determined independently of the third Ke- 
pler law: P, K, q from optical spectroscopy of HDE 226868, 
and m from the analysis of certain properties of the Cygnus 
X-1 X-ray emission. 
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OA// sin i, where i is the inchnation angle of the orbital 
plane to the plane of the sky, as 
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An import a nt res ult of the spectroscopical analysis 
bv lGies et al.l (|2003l ) is that the emission component of 
the radial velocity is compatible with a mass ratio 



q= — =0.36 ±0.05. 

M 



(2) 



The recent result by Shaposhnikov fc Titarchuk ( 2007t ) 
for m yields for the mass of HDE 226868 



M= (24±5)Mo; 



(3) 



it disagrees with M = 17. SM© bv lHerrero et al. 
which, among other things, relies upon various es- 
timates of i. Even larger is the disagreeme n t with 



Ziolkowskil |2005l ) 



M = (40 ± 5)Mo obtained by 
with evolutionary calculation for the most likely in- 
tervals of the values of the distance and of the ef- 
fective temperatur^l]: 1.95 to 2.35 kpc and 30000 to 
31 000 K. It is important to note t hat the method 
by IShaposhnikov fc TitarchukI (|2007l ) is not only in- 
dependent of the lingering uncertainty in the system's 
distance from us, but it was also robustly tested by 
reproducing the masses of other black hole candi- 
dates, hke_the_one_Jn_them GRS 1915 ± 
105 (ICastro-Tiradq et al.l 1 19921 ) , obtained with IR 



( Greiner et al 



200^ and X-ray jshrader fc Titarchuk 



20031) observations. Thus, in the following we will trust- 
worthy rely upon such a result. Note that the mass 



estimate by IShaposhnikov fc TitarchukI (|2007l ) for m is 
we ll within the range 4 .8Mo < m < 14.7M0 obtained 
by Herrero et al. ( 19951 ) . 



We are now in the position of looking for genuine 
deviations from the third Kepler law by comparing the 
phenomcnologically determined period P to the calcu- 
lated Kcplerian one which, from eq. ([1]) and 



XM 



q J sm I 



(4) 



for the relative semimajor axis, can be expressed as 



pKop 




(5) 



in terms of quantities determined independently of the 
third Kepler law itself: note that eq. (O does not 
depend on M. The overall uncertainty 5P^'^^ in eq. 
([5]), which will be considered as a function of z, can be 
conservatively evaluated by linearly adding the various 
mismodelled terms due to (Jq, (JG, Jm, (5 A', 8P as 
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In Figure |T] we plot the ratio 

pKcpl 



AP 



IP 



5(AP) SP + 5PKcp 



(7) 



as a function of the inclination; outside the range 41 
deg ^ j ^ 56 deg, in which the purely Keplerian period 
is an adequate model for the measured orbital one, de- 
viations from the third Kepler law significant at more 
than 1-sigma level would occ ur; in part i cular, the es- 



timate i = (30 ± 7) deg by iGies et al.l (|2003l ) would 
yield a discrepancy signifi cant at 4/2-sigma-lev el, while 
for z = 35 deg, chosen by Herrero et al. ( 19951 ) because 
lying in the middle of almost all determinations in lit- 
erature, the deviation from the third Kepler law would 
amount to 2-sigma. Can such deviations be considered 
physically meaningful? 

In principle, also the first post-N ewtonian (IPN) cor- 
rectio n to the Kcplerian period ( Damour fc Deruelld 



se pa 

and 28000 to 32000 K. IZiolkowskil 1^0^) obtained a range 29-^-50 
for M. 



P 



IPN 



V2c2 



1 



nGmKP 



3(i + g)2 



(8) 




Fig. 1 Discrepancy between tlie plienomenologically mea- 
sured orbital period P and the computed Keplerian one 
P^'^'P as function of the inclination i. 



here expressed in terms of the phenomenologically de- 
termined quantities, should be taken into account in 
modeling the orbital period; however, it is not possible 
that eq. ^ is the cause of the significant discrepan- 
cies because it is six orders of magnitude smaller than 
jpKcp j^ggif fgj. oil values of i. Note that considering the 
possibility that IPN terms may have some observable 
effects in so close binary systems should not be consid- 
ered trivial since for, e.g., the exoplanet HD209458b the 
accuracy in determining its 3.5 d orbital period would, 
in fac t, allow to detedljthe relativistic correction I loriol 
20061) . 

Let us see what happens if we add to the Keple- 
rian period the correction due to the quadrupole mass 
moment Q of the system: in terms of the determined 
parameters we have, from eq. ()23|) . 



^_ k2 ll'mPyCK JUK^ I KP 



3\/2 V V / ^'1^* 6 ^TTsini^ 
where we have posed 

_ 1 + g 



(9) 



(10) 



^The uncertainties in the system's parameters induce systematic 
errors which prevent from implementing such a goal. 



Fig. 2 Upper panel: Q/5Q. Lower panel: Q/MB? as 
functions of the inclination i. 



The total uncertainty in Q can be conservatively eval- 
uated by summing the various contributions 



5Q\g < 



5Qln < 



5Q\p < 



5Q\k< 
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(11) 



In Figure [Hwc plot \Q\/5Q (upper panel) smc^Q/MH^ 
(lower panel) : we see that ascribing the deviation from 
the third Kepler law, when it becomes significant, to 
the HDE 226868 's quadrupole mass moment leads to 
meaningless results since Q would assume values in- 
compatible with zero at more than 1-sigma level, but 
unreasonably large. 

By conservatively assuming that no deviations from 
the third Kepler law occur, given the present-day level 
of accuracy in knowing the system's parameters, it is 
possible to determine i as 



arcsin KP^^^ ( 2T:GmK 



-1/3- 



= (48 ± 7) deg. 



^We used R = IJRq jHerrero et al.l Il995h : according to 
IZiolkowskil bOO^ . R/Rq = 10. SM, with d in kpc: doil.S- 2.35 
kpc. 
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Fig. 3 Relative semimajor axis a, in Solar radii, as function 
of the inclination i. 



Fig. 4 Mass of HDE 226868, in Solar masses, as function 
of the inclination i. The third Kepler law was used. 



(12) 



In regard to the inclination i, Gies et al. ( 20031 ) ob- 
tain 



i = (30 ± 7) deg 



(13) 



from the measured emission and absorption radial 
velocities and by making certain assumptions about 
isotropy of the outflow velocity; after noting that their 
result is compatible with other estimate^ like 28 deg 
<: i < 38 deg from light curve analysis bv ICies et al. 
|l986al) and 10 deg < i < 40 de g from fitting the 



X-ray emission curve bv IWen et al.l (|l999r ). iGies et al. 
I 2OO3I ) conclude that eq. (|13p is likely a lower limit 



Our result , which is incompatible with the one by 



e.g., the result bv iDolan fc Tapial (|l989r ). can be use- 
ful in terms of anisotro py in the wind outflow velocity 
( Friend fc Castoilll982h . 

It is also interesting to evaluate the relative sep- 
aration of HDE 226868 and Cygnus X-1 (see Figure 
[3]): such an information may be helpful about HDE 
226868's dimensions. The relative semimajor axis, ac- 
cording to cq. ([T]) and cq. ^ which are independent of 
the third Kepler law and on the mass m of Cygnus X-1, 
amounts to a (42 ± 9)Rq (w 0.2 AU) for i ^ 48 deg; 
the minimum value, i.e. a — 31.6i?0 (0.1 AU), would 
occur for the edge-on configuration while for i = 10 deg 
it amounts to about 181 Solar radii (0.8 AU). 

Note that, should one consider M = (40 ± 5)Mq 



inf erred bv IZiolkowskil (|2005|) . the black hole's mass 
by IShaposhnikov fc TitarchukI (j2007l ) would yield q = 



^Optical polarimetry of HDE 226868 yields 25 d eg < i < 67 d 
jPolan fc Tapialll989h and 20 deg < i < 40 deg llDanie]||l98lh . 



0.22 ± 0.05. Such a value, w hich is inconsisten t at 1.4- 
sigma level with the one by Gies et al. ( 2003[ l. would 



lead to unreasonably large values for the quadrupole 
mass moment for « ^ 45 deg, being all the other values 
up to 90 deg compatible with zero Q and no devia- 
tions from the third Kepler law. Near edge-on orbital 
configurations imply that eclipses in the X-ray emis- 
sion of Cygnus X-1 would be allowed {i > 65 deg), but 
they have never been detec ted. Thus, t he es timate for 
the HDE 226868's mass bv IZiotkowski (|2005l ) does not 
seem plausible. Another line of reasoning is the fol- 
lowing. From the third Kepler law we can express the 
HDE 226868's mass in terms of the spectroscopically 
determined parameters as 



Gies et al.l ( 20031 ), being, instead , com patible with, ^^'^ 



qsim 



(14) 



as can be noted from Figure HI M assumes the values 
bv IZiotkowskil I2OO5I ). who relies on the third Kepler 
law, for 37 deg ^ * ^ 41 deg; but our previous analysis, 
independent of M itself, has shown that for such values 
significant deviations from the third Kepler law would 
occur (see Figure [T]). 

Note that our results ar e compatible wit h a system 
distance of about 1.80 kpc (lMalvshevalll997l) and an ef - 



fective temperature Tg — 32000 K (jHerrero et al.lll995l) 
of HDE 226868, according to Table 1 of dZiolkowski 



2OO5I) jased on the quite general approach by lPaczynskil 
( I974I) . 



3 Discussion and conclusions 

In this paper we refined some orbital and physical pa- 
rameters of the HDE 226868/Cygnus X-1 binary sys- 
tem by exploiting the latest spectroscopic optical data 
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from HDE 226868 and the recent determination of the 
mass m of Cygnus X-1 via QPO-index correlation scal- 
ing. For the mass of the supergiant star we obtained 
M = (24 ± 5)Mq, which is in disagreement with other 
estimates present in literature. Then, we constrained 
the inclination i by looking for statistically significant 
deviations from the third Kepler law which would oc- 
cur, at more than 1-sigma level, for i < 41 deg and 
z > 56 deg. After noting that they cannot be due to the 
IPN correction to the Keplerian period, we concluded 
that, if attributed to the star's quadrupole mass mo- 
ment Q, they are to be considered unphysical because 
of the too large values which Q would take in this case. 
The inspection of the possible values of Q allowed also 
to consider unhkely the value M = (40±5)Mq for HDE 
226868 recently inferred from evolutionary calculation 
because it would allow for high values of i compatible 
with eclipses which, instead, are so far absent. Thus, 
by reasonably assuming that the Keplerian period is 
an adequate model of the orbital period we determined 
i = (48 ± 7) deg. In this case, the relative semimajor 
axis amounts to a = (42 ± 9)Rq (w 0.2 AU). 



Appendix: the correction to the orbital period 
due to the quadrupole mass moment 

One of the six Keplerian orbital elements in terms of 
which it is possible to parameterize the orbital motion 
in a binary system is the mean anomaly A4 defined as 
M = n{t — Tq), where n is the mean motion and Tq 
is the time of pericenter passage. The mean motion 
n = 27r/Pb is inversely proportional to the time elapsed 
between two consecutive crossings of the pericenter, i.e. 
the anomalistic period Ph- In Newtonian mechanics, 
for two point-like bodies, n reduces to the usual Ke- 
plerian expression n^'^P = a/ GM/a^, where a is the 
semi-major axis of the relative orbit and M = mi -I- m2 
is the sum of the masses. In many binary systems the 
period Ph is accurately determined in a phenomcno- 
logical, model-independent way, so that it accounts for 
all the dynamical features of the system, not only those 
coming from the Newtonian point- like terms, within the 
measurement precision. 

Here we wish to calculate the contribution of the 
quadrupole mass moment Q to the orbital period in a 
general way. 

By assuming that the spins and the orbital angu- 



lar m omentum are aligned, as done in (jHerrero et al 



19951) for the HDE 226868/Cygnus X-1 system, the ra- 
dial component of the acceleration induced by Q can 
be cast into the form 



Ac 



3GQ 



while the other component, i.e. the latitudinal one, van- 
ishes. Aq is small with respect to the usual Newtonian 
monopole, so that it can be treated perturbatively. In 
order to derive its impact on the orbital period Pb, let 
us consider the Gauss equation for the variation of the 
mean anomaly in the case of an entirely radial disturb- 
ing acceleration A 



dM 2 ^ 

— — =71 A 

at na 



(1 



-A cos /, 



(16) 



where / is the true anomaly, reckoned from the peri- 
center. After inserting Aq into the right-hand-side of 
eq. (|16p . it must be evaluated onto the unperturbed 
Keplerian ellipse 



1 + e cos / ' 



By using (|Rovll2005D 



eq. pB)) yields 



(17) 



(18) 



dt 



n{l + ecos/)^ 
(l-e2)3/2 



3Gg(H-ecos/)3 

2\3 



COS /(I -I- ecos/) 



2e 



n'^w^{l — e^) 



- 1 



(19) 



The orbital period can be obtained as 



Ph 



(l_e2)3/2 



2tt 



1 - 



cos /(I -I- ecos/) 



3Gg(l-Fecos/)3 
n^a^(l — e^)3 
df 



2e 

From eq. (|20p it follows 
Ph = P^°^ + P'^, 
with 

pKcp _ o / 



1 



(1 + e cos /)■ 



r(20) 



pQ ^ 



(21) 



(22) 



^GaM3(l-e2)3 



Note that P'^ in eq. ([22|) agrees with the expression of 
the anomalistic pe riod of a satellite orbiting an oblate 
planet obtained by Capderoul ( 2005^ ): for a direct com- 
parison Q = — M_R^J2; where J2 is the first even zonal 
harmonic of the multipolar expansion of the Newtonian 
part of the gravitational potential of the central body 
of mass M and equatorial radius R. 



(15) 
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Solving for Q, wc finally get 
Q^^V ^ L_ Ma2(l-e2)3/2. (23) 

OTT O 
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